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Abstract 

In a typical exclusive quarkonium production process, when the center-of-mass energy, ^/s, is 
\ much greater than the heavy quark mass m, large kinematic logarithms of s/m^ will unavoidably 

arise at each order of perturbative expansion in the short-distance coefficients of the nonrelativistic 
QCD (NRQCD) factorization formalism, which may potentially harm the perturbative expansion. 
This symptom reflects that the hard regime in NRQCD factorization is too coarse and should be 
further factorized. We suggest that this regime can be further separated into "hard" and "collinear" 
^ ■ degrees of freedom, so that the familiar light-cone approach can be employed to reproduce the 

NRQCD matching coefficients at the zeroth order of m^/s and order by order in ag. Taking two 
simple processes, exclusive ijb + 'j production in e^e~ annihilation and Higgs boson radiative decay 
into T, as examples, we illustrate how the leading logarithms of s/m^ in the NRQCD matching 
coefficients are identified and summed to all orders in with the aid of Brodsky-Lepage evolution 
equation. 
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Introduction. One of the classical applications of perturbative Quantum Choromo dynam- 
ics (QCD) is the successful description of many exclusive processes with large momentum 
transfer using coUinear factorization, that allows one to express the scattering amplitude as 
the convolution of perturbatively calculable short-distance parts and the nonperturbative 
but universal light-cone distribution amplitudes (LCDA) jl|, Q] • In describing hard exclusive 
processes involving heavy quarkonium, i.e., a nonrelativistic bound state made of a heavy 
quark and a heavy antiquark, however, there also exists another widely-accepted theoretical 
framework, the nonrelativistic QCD (NRQCD) factorization formalism j^. In this approach, 
the amplitude can also be put in a factorized form, that is, an infinite sum of products of 
short- distance coefficients and nonperturbative, albeit universal, NRQCD matrix elements. 

In recent years, considerable amount of efforts have been spent to understand exclu- 
sive charmonium production mechanisms from both NRQCD and light-cone perspectives. 
This endeavor is perhaps largely propelled by a somewhat unexpected finding, that the 
lowest-order NRQCD calculation of double charmonium production rate for the process 
e+e^ J/iprjc fell short by about one order-of-magnitude of the BELLE measurement j^. 
The validity of applying NRQCD to charmonia was soon questioned by some authors, who 
advocated that if charmonium is treated as light mesori, the light-cone approach instead 
could satisfactorily accommodate the Belle data 0,0,0]. However, a careful reexamina- 
tion 9j suspected that these optimistic assertions are premature and it was argued that a 
"consistent" light-cone analysis would in fact yield a result not much different from NRQCD, 
so that the situation has not truly improved. 

In reviewing this episode, one may get the impression that the light-cone and NRQCD ap- 
proaches are two drastically different, and, competing, theoretical frameworks. Indeed, these 
two approaches are rooted in two different types of operator-product-expansions (OPE), 
where the former is intimately linked to the light-cone (twist) expansion, and the latter is 
more closely related to a local (large mass) expansion. 

The main motif of this work is to show that, these two approaches, both bearing solid 
theoretical grounds, need not to be regarded solely as rivals. As a matter of fact, they can 
benefit each other ^. In particular, it turns out that the concepts and techniques developed 



^ We note that there already exist some attempts along this line, but motivated by somewhat different 
considerations. In particular, there is a recent work trying to bridge the leading-twist LCDA of S'-wave 
charmonia with the local NRQCD matrix elements [l^ (see also [HI. 
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in the light-cone approach can be fruitfully exploited to improve the NRQCD calculations. 
For example, we will argue that, for some class of exclusive single- quarkonium production 
processes, the NRQCD short-distance coefficients can be reproduced order by order in a^, in 
the language of collinear factorization. Moreover, the light-cone approach can easily resum 
large kinematic logarithms appearing in the NRQCD short- distance coefficients. 

Refactorization of NRQCD matching coefficients. NRQCD is an effective field theory 
that is constructed to portray the nonrelativistic character of heavy quarkonium. The slow 
relative motion between quark and antiquark in a quarkonium, v, naturally constitutes the 
expansion parameter of the theory. To produce a quarkonium, the heavy quark pair must 
be created within a distance shorter than 1/m and have a small relative velocity, to warrant 
a significant probability to form a quarkonium. The first condition guarantees that the 
parton process can be calculated perturbatively owing to asymptotic freedom. The latter 
requirement implies that the full amplitude can be expanded in power of v, and NRQCD 
endows a well-defined meaning for this Taylor expansion procedure. 

For definiteness of our discussion, we will confine ourselves to the bottomonium throughout 
this work, whenever a quarkonium is referred to. Of course, all the results can be trivially 
copied to the charmonium case. Furthermore, for simplicity, in this work we will concentrate 
on the single quarkonium production process, and not touch upon the double-quarkonium 
production for its additional technical complication. For a typical hard process that creates a 
bottomonium, generically denoted by H, NRQCD factorization demands that the amplitude 
can be expressed as 

M[H] ~ J]C„(Q,m)(i7|a|0), (1) 

n 

where On stands for an appropriate local color-singlet NRQCD operator, m is the b quark 
mass, Q refers to a typical kinematic scale such as the center-of-mass energy, presumably 
much greater than m. The important feature is that each On has a definite power counting 
in V, so that the expansion in ([1]) becomes practically maneuverable. 

NRQCD factorization manifestly separates the effects of hard degree of freedom (p^ > m) 
from the remaining lower-energy ones ^. In Eq. ([1]), the matching coefficients Cn encode the 



^ There are totally three kinds of such, soft (p'' ~ mv) , potential (p^ ^ mv'^, p ^ mv) and uUrasoft 



(p'^ ~ mw^) 13|. All these scalings are specified in the quarkonium rest frame. 
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effects of hard quantum ffuctuation, which are perturbatively calculable and infrared finite, 
whereas the local NRQCD matrix elements characterize the long-distance, and primarily 
nonperturbative, aspects of H. 

The short-distance coefficient C„ depends on two energy scales, Q and m. When these 
two scales are widely separated, it can be naturally organized by Taylor series in m^/Q^, 
practically only the leading power term needs to be retained ^. It is always possible, by 
factoring out an appropriate kinematic factor proportional to some powers of Q, to render 
Cn dimensionless, and such that it only depends on th? and through their dimensionless 
ratio and it scales with m?'/Q'^ at most logarithmically. 

It is conceivable that at each order in loop expansion, C„ receives contribution from 
logarithms of m^/Q"^ as well as constants. When Q ^ m, there arises the ambiguity of 
setting the renormalization scale in the as, whether to put it around or may result 
in quite different answer for a fixed order calculation. Moreover, these large logarithms 
may potentially ruin the perturbative expansion in strong coupling constant. Solution to 
these problems lies beyond the scope of NRQCD. Needless to say, it is highly desirable that 
these large kinematic logarithms in the NRQCD matching coefficients can be identified and 
summed to all orders in a^, in order to ameliorate perturbative expansion. 

The cause for this symptom can be easily traced. The problem is that the hard quanta 
integrated out by NRQCD still contains two widely separated scales, accommodating the 
quantum fluctuations extending from the scale m to a much higher scale Q. It is natural to 
speculate, for such a multi-scale problem, whether one can divide this rather coarse region 
further into finer ones. To achieve this, one needs first identify various relevant degrees of 
freedom, then disentangle their contributions, and finally express C„ in a factorized form. 

The NRQCD factorization is justified mainly because b is heavy, m ^ Aqcd- As Q ^ m, 
the H moves almost with the speed of light in a natural reference frame such as the center- 
of-mass frame. Since the heavy b and b move nearly along the light cone, they can be 
considered to be light in the kinematic sense. We thus are facing a challenge where the 
conflicting nature of "light" heavy quarks must be consistently incorporated. 

From now on, we will specialize to the S-wave quarkonium production. For simplicity. 



^ In fact, for some classes of exclusive charmonium production processes, it has been proved that NRQCD 
factorization may no longer hold beyond the leading order in /Q^ lj| . Therefore it seems not profitable, 



for both pragmatic and theoretical reason, to proceed beyond the leading power. 
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we will only consider the leading term in the NRQCD expansion, since Oi is the simplest 
NRQCD operator and constitutes the most important contribution. However, we would like 
to stress that the purpose of imposing this restriction is mainly for simplicity. Since our 
reasoning will be based on quite general ground, there should be no principal difficulty to 
proceed to higher orders in v expansion. 

To compute the short- distance coefficient Ci, a convenient method is to consider the quark 
amplitude where H is replaced by a free bb pair sharing the same quantum numbers (color, 
spin, orbital angular momentum) as the leading Fock component of H. At the lowest-order in 
V, b and b are forced to partition equally the total momentum of the pair, P. For definiteness 
P is supposed to move fast along the positive z axis. Consider a generic loop diagram that 
contributes to the corresponding quark amplitude. The so-called "method of region" 15 1 
can be utilized to help identify the relevant degrees of freedom in Ci. There is a ^^hard" 
region in which the loop momentum scales as ~ Q. In this region, the b quark can 
be treated as massless. There are several kinds of infrared modes, soft {jf ~ m), collinear 
(p+ ^ Q,p^ ~ m?/Q,p±^ ~ m), and anti-collinear (p+ ~ vn?jQ^p~ ^ Q,p± ^ m) ^. All 
these infrared modes have the virtuality of order m^, hence the mass of b quark must be 
retained in these infrared regions. The validity of NRQCD factorization guarantees there is 
no any overlap between these "infrared" quanta and those truly infrared modes intrinsic to 
NRQCD such as the potential mode, as can be readily distinguished by their virtualities. 

It turns out that at the leading power of m^/Q^, the soft and anti-collinear modes do 
not contribute to Ci p^. Therefore the relevant dynamic degrees of freedom are only the 
hard and collinear modes, which resembles the hard exclusive processes involving a single 
hadron such as 71-7 transition form factor [1]. Quite naturally, one can follow the standard 
collinear factorization theorem at leading twist to cast Ci as 



where T refers to the hard-scattering part involving massless quarks, and can be viewed 
as the LCDA of a free bb pair with the same quantum numbers as H (Later we will give 

^ We caution that the term hard here only refers to a subset of hard region in the sense of NRQCD, and 
the actual meaning of this term may vary depending on different context. Similarly the term "soft" may 
also bear different meaning in different places. 

^ Our convention in defining the light-cone variable is — ^ ^ , and a four- vector is decomposed 
into A^) in this convention. 
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a rigorous definition for tliis LCDA). x is tlie fraction of plus- momentum carried by the b 
quark with respect to that of the pair, and the operation ® imphes the convolution over x 
between two parts, fi is an arbitrary scale separating these two regions. It is important to 
note that the scales Q and m now become fully disentangled, i.e. the T only depends on Q 
but not on m, likewise the cannot have any direct sensitivity on the hard scale Q. 

It is worth emphasizing that introduced here (to be concrete, at this moment let us 
specialize to the twist-2 LCDA of the bb{^S^\P) state), bears some important difference 
with the twist-2 LCDA of a pion. The latter is a genuinely nonperturbative object that is 
sensitive to the quantum fluctuation at hadronic scale, so that it can only be extracted from 
experiments or from some nonperturbative tools such as lattice simulation or QCD sum 
rules. On the contrary, m acts as the infrared cutoff in the former case. Since m is much 
greater than Aqcd, thus can be reliably computed in perturbation theory. However, we 
note that, since both LCDAs of the bb{^Si^^\ P) state and a pion possess the same ultraviolet 
behavior, they obey the same renormalization group (RG) equation. 

The factorization formula ([2]) is expected to hold to all orders in as- Both of the two 
ingredients, T and are perturbatively improvable. When Q ^ m, this light-cone-based 
method provides an alternative and efficient means to reproduce the NRQCD short- distance 
coefficient, including the logarithms as well as the constants, to any desired order in as- 
At tree level, the hard part T{x, Q) can be easily evaluated, and the defined at the scale 
~ m is simply 5{x — |), to be compatible with the NRQCD matching condition. Using 
the standard evolution equation to evolve this LCDA into higher scale /i ~ Q, combined 
with the tree-level result of T, Eq. ([2]) then automatically sums the leading logarithms to 
all orders in ag- 
in the remainder of the paper, we will test our understanding through two explicit ex- 
amples of exclusive production of a bottomonium plus a photon. We will illustrate how the 
leading kinematic logarithms in the NRQCD short-distance coefficient can be identified and 
resummed within the framework of the factorization formula ([2]). 

Example 1 — e"'"e~ rjb'y. Let us first consider the exclusive rjb + ■y production at high- 
energy electron-positron collision, say, at LEP experiment. Similar processes at B factory 
environment have been studied in a light-front model 17| and the NRQCD approach To 
our purpose, it suffices to focus on the decay of a highly virtual photon, YiQ) ~^ Vb{P)+l{k). 
According to the NRQCD factorization, one can express the decay amplitude at the lowest 



6 



order in v as 



M[Y ^ + 7l = M[Y B(\Sj,'\P) + %^ , (3) 

where ip'^ and x ^i-re two-component Pauli spinor fields. The factor associated with the 
NRQCD matrix element is chosen for convenience, mainly to compensate the fact that the 
particle states in the amplitude is relativistically normalized, but the r?^ state in NRQCD 
matrix element conventionally assumes non-relativistic normalization [S]. It is worth men- 
tioning that, the factor {rji^ip'^ x\^) / coincides with f^^, the r^^ decay constant, to the 
lowest-order accuracy in as and v. 

M.^ being the short-distance coefficient, represents the quark amplitude when rji, is re- 
placed by the free h and h quarks that equally share the total momentum P, and have the 
spin-color wave function l^^-^-^^-*-^^ x where A^^^ = 3 is the number of color and j repre- 



sent color indices. The calculation is usually expedited 
method 11911, though alternative methods also exist 120 



utilizing the covariant projection 



2l|. The result is 



19|, though alternative methods also exist 

A?[7* - +7] = v^^6,.„^4e7Q"fc^F , (4) 

where = — | is the fractional electric charge of h quark. We introduce F as a dimensionless 
form factor. It can be simultaneously expanded in powers of and rn? jQ^: 

where Cni are constants. As explained before, the "higher twist" contributions suppressed 
by l/Q^ are omitted. For this simple process, the leading logarithm at the nth loop order 
is of the form a" ln"'((5^/m^), and our goal is to deduce C„„ analytically as well as resum 
their effects to all orders in a<( ®. 

The difficulty of NRQCD matching computation rapidly grows once beyond the tree 
level, and the resulting analytic expressions of F are generally rather cumbersome, often 
plagued with special functions depending on m? /Q"^ in a complicated way (for some concrete 



examples on single and double charmonium production to one-loop order, see Ref. 



22 



^ Note that the function F is complex- valued in this process. Its imaginary part can contain sub- leading 
logarithms. However, insofar as the leading logarithms are concerned, we are allowed to ignore the 
contributions from the imaginary part. 
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24l|). However, the expanded form as indicated in ([5]) is always very simple. It is desirable if 
F can be reproduced in a more efficient way. As discussed earlier, this is indeed possible, if 
we try to disentangle "hard" and "collinear" quanta in the NRQCD matching coefficients. 
For this example, a factorization formula in complete analogy to the tt— 7 form factor is 
expected to hold: 



(6) 



where /i is an arbitrary factorization scale that separates the hard and collinear contributions. 
The hard function T is identical to what appears in the 71—7 form factor, which can be 
expanded as T = T^o) + ^^T^^) + ■ ■ ■ . The tree-level result^(°l(a;, Q) = ^ + ^z^, is well 
known, and T*-^-* has also been available for a long time 25|, l26|]. The "long-distance" 
contribution is characterized by 0, the twist-2 LCDA of a free hb pair in the spin-color state 
^5*0^^ and with total momentum P halved by its two constituents. This LCDA admits a 
definition in term of operator matrix element 2|: 

1 f dw' 



m, yu) 



-ixP~^i 



i4^\P)|6(0,^^,0^)7+75&(0)|0), (7) 



V2NcJ 27r 

where for simplicity we have suppressed the light-like gauge link. This definition is boost 
invariant along the z axis. The can be systematically expanded perturbatively, = 
m, /i) + m, /i) + ■ ■ ■ , here /i can be identified with the renormalization scale 

of the operator matrix element. Its tree-level result is very simple, (p^^^x, /i ~ m) x) = 

6{x — ^), embodying the NRQCD matching condition. For the investigations on (f)^^\ one 



may consult Ref. 10 



Substituting the tree-level expressions of T(°) and 0^°) into 



one then obtains Cqq = 4, which agrees with the result derived from the standard NRQCD 



matching 18j. 

The invariance of physical amplitude about the choice of the factorization scale /i leads to 
the RG equation, which is conventionally referred to as Brodsky-Lepage (BL) equation 2?! . 
This equation governs the evolution of the LCDA from one scale to another one, through 
which the collinear logarithms can be resummed. The BL equation reads: 



_d_ 



■0(x,/i) 



TT 



c 

dy^Vo{x,y)(i){y,fi), 



where 



Vo{x,y) 



X 



i-y 



1 + ]e{x-y) + - 

x-yj y 



1 + 



y-x 



e{y 



X) 



(8) 



(9) 



is the kernel governing the evolution of LCDA of a helicity-zero meson. Cp = is 
the Casmir for SU (Nc) fundamental representation, the subscript "+" implies the familiar 
"plus" prescription. This kernel can be interpreted as the amplitude for a quark with 
fractional plus-momentum x and antiquark with 1 — x to become a quark with fractional 
plus-momentum y and antiquark with 1 — y hj the exchange of one collinear gluon. 

As is well known, the kernel Vq admits the eigenfunctions G'„(x) = x(l — x)Cn^'^^ {2x — 1), 
which are Gegenbauer polynomials of order | multiplied by the weight function x(l — x): 

I 

dyVo{x,y)Gn{y) = -7„G'„(x), (10a) 



n+l 

i 

In 



1^1 1 . , . 

- + 2> -. 10b 

2 (n+l)(n + 2) ^ ' 



It is a standard procedure to decompose the 0(x) in the basis of Gn- 

oo 

0(x,/i) = ^0„(/i)G„(a;), (11a) 



n=0 

•1 



0„(^) = 4(2n + 3) , dxC^^/^\2x-l)^{x,fi). (lib) 



(^ + l)(^ + 2) Jo 

At leading logarithmic accuracy, one can express the form factor as the infinite sum of 
Gegenbauer moments: 

^ = / rfxr(°)(x)0(°)(x,g^) = 5:ff(Q^). (12) 

All the odd Gegenbauer moments vanish. In deriving f|T2l) . we have made use of the fact 
dx C2^^'^\2x — 1) = 1 for any integer n. Using C2j'^\o) = ^—^^j^^-^^ and evolving the 
the moments from the low scale /i ~ m to Q, we obtain 

'(5).=§*-(3a)*-. 

4(-l)"(4n + 3) ):'^ (13b) 



= 4(-l)"(4n + 3)i^^' 
2" ^ ^ ^ ^(2n + 2) 



The anomalous dimension for each moment is dn = 2CFln/ l^o-, where 7„ is given in fllObp . 
/3o = yA'c ~ l'^/ is the one-loop coefficient of QCD [3 function and ^/ is number of active 
quark flavors. The running strong coupling constant is given by as{Q'^) = in(Q2/A^ — y. (t>'p 
is an alternating series and its magnitude declines very slowly (oc 1/ y/n asymptotically 11|). 
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and the convergence of the sum becomes extremely slow. Practically, one has to include an 
extraordinary number of terms in the sum to warrant the convergence of the result. 

In Fig. [T]we show the effect of leading logarithm resummation, fll3al) . for a wide range 
of values of Q. One can clearly see that this RG-improved result seems to only have a 
minor numerical impact, perhaps less important than the constant terms in the fixed order 
calculation. As oo, only the first moment in fll3al) can survive since c/q = and all 

the other d2n are positive, therefore F — > |Coo- However, as might be inferred from Fig. [U 
to reach this asymptotic result, an unphysically large Q {i.e. ^ Mpianck) is required. 

Even though leading logarithms have been summed to all orders in as, it would still 
be illuminating if the coefficient of each leading logarithm at n-th loop order, C„„, can be 
deduced. To this end, it is more convenient to work in the x space than the moment space. 
We can rewrite F as: 

fI"^] = [\xT^^\x)^'^''\x,Q)= [\xT^^\x) expire fVo*]^'^''\x), (14) 

\Q J LL Jo Jo 

where 

In the second equation of Eq. f|T^ . we have substituted the symbolic solution of the BL 
equation for (j){x,Q), which is evolved from the infrared scale at ;U ~ m to the high scale 
fi ^ Q. The operation implies that, upon the expansion of the exponential, the result- 
ing products of the argument of the exponential are to be treated as convolutions. It is 
worth noting that, Eq. ( |T4|) is essentially the statement that the moment-space amplitude 
exponentiates, as indicated in fll3ap . 

We can expand equation (IMI) iteratively: 

F ("^^ = f'dx (x) (x) + kCf f'dxf'dy r(o) (x) Vo{x, y)^^^^ (y) 

\Q / LL Jo Jo Jo 

+ ^ fdxfdyfdzT^'\x)V,{x,y)V,{y,z)^^'\z) + .... (16) 

^! Jo ^0 Jo 

Equation ( |T6i) admits a clear physical picture how the leading logarithms are built up. At 
the hard vertex, h and h are produced to carry arbitrary fractional plus-momentum x and 
1 — X, and transverse momenta of order Q. In the course of evolution, they keep exchanging 
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collinear gluons ^ and reshuffle the respective fractional plus-momenta, and gradually de- 
crease the transverse momenta. When ready to transition to a physical rjb state, they finally 
adjust themselves to become nearly collinear and halve the total momentum. 

We choose the order of the multiple integration in flTB]) from the left to the right, and 
leaves the integration over the (j)^'^^ in the last step. With the help of (lAlaP and (lAlbl) . we 
are able to deduce the leading logarithms through the two-loop order: 



00 



LL 



+ Ci 



1 + 

(47r)2 



Ait 
InM'^ 



(17) 



All 



In 2 



81n2 + ln^2 



+ 



where the coefflcient Cu agrees with the asymptotic expansion of the exact NLO correction 
to the NRQCD short-distance coefflcient |28|], and can also be read off by imposing x = ^ 
m T(i)(x) in Ref. [26|. Our prediction to C22 is new. From Fig. [1], one can see that, for a 
wide range of Q/m, the leading logarithm approximation through two- loop order is already 
quite close to the resummed results, ( 113ap . 
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FIG. 1: The form factor F{rin? /Q'^) as a function of Q for fixed mj, = 4.7 GeV. The left panel 
is for 7* — > ryj, + 7, and the right panel for /i — > T + 7. We have fixed ra/=5, /3o = ^ and 
^QCD = 100 MeV. Various curves represent the improved results to different extent in leading 
logarithm approximation. The dot-dashed line only includes the logarithm to one- loop order, the 
dashed line includes the leading logarithms through two- loop order, and the solid line sums the 
leading logarithms to all orders in a^- 



^ We have tacitly made an important simplification. The kernel Vb we use is valid only for the flavor 
nonsinglet where annihilation of qq is not allowed. In our case, it is possible that the quark pair in the 
intermediate loop to annihilate into two gluons, then recreates hh pair. These types of singlet diagrams 
first occur at two loop. In this work, we neglect the logarithms arising from such singlet contributions. 
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Example 2 — Higgs radiative decay into T. It is tempting to test our understanding in 
a slightly different situation, e.g. a process that a transversely polarized T is produced at 
leading power. Let us consider h{Q) — > T(P) +7(A;). In view of its clean signature, this 
decay channel may be potentially useful to search for the Higgs boson in the forthcoming 
CERN Large Hadron Collider (LHC) experiment, were it not completely swallowed by the 
much more copious background. In accordance with the NRQCD factorization, we can 
express the amplitude at the lowest order in t> as ^ 

A1|A ^ T + ,1 = Mil, U.CS^Kp) + <^W^;;-^^I°>. (18) 

It might be worth pointing out that, the NRQCD matrix element {T{e)\ilj''(T ■ ex\0)/y/rn 
coincides with /x, T decay constant, to the lowest-order accuracy in as and v. 

Analogous to the first example, Ai here represents the parton amplitude when T is 
replaced by the free b and b quarks in the ^S[^^ state, each of which carries half of the total 
momentum P. It can be expressed as follows: 

M[h^bbCsi'\P)+^] = -Vm'-Y9H,ie*r-e;F(^^y (19) 

where = M|, gf^^f^i = m/v is the Yukawa coupling between Higgs and b quarks, and 
V = 246 GeV implies the vacuum expectation value of the Higgs field. Since the photon 
must be transversely polarized, so is the T. 

The dimensionless form factor F here is expected to also satisfy a collinear factorization 
formula similar to (Q. The tree- level hard coefficient function T^^-* is the same as the 
preceding example, but here one should use the twist-2 LCDA of the transversely polarized 
bb{^S[^\P) state, which can be defined as the following operator matrix element [2]: 

' V2KJ 271 

X {bbCs['\P, A^)| 6(0, W-, 0^) a+%(A^) 6(0) |0) . (20) 

At the lowest order in a^, (j)f^{x.,jji ^ m) = (f)^^\x) = S{x — 1/2). Thus in this example, we 
also get Coo = 4. The (f)± also satisfies the BL equation ([S]), except the kernel Vq should be 



Our main concern here is for the illustrative purpose, so we do not target at a serious phcnomcnological 
analysis of this decay channel. For instance, we have not included the fragmentation contribution from 
h 7*7 — > T7. 
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replaced by a different one: 



V±{x,y) = Vo{x,y) - 



^ "'x-y) + -9{y - x) 



(21) 



.1-2/ y 

This new Icernel permits the same eigenfunctions Gn{x) as given in Eq. fllOap . but the 
corresponding eigenvalues are slightly different from Eq. fllObI) . i.e. 7„ = i + 2 X]j=2 j- 

Completely analogous to the preceding example, we can express the form factor F as the 
infinite sum of Gegenbauer moments, with evolution effect taken into account: 

(ia;T(o)(x)0f (x,Q') 



LL JO 



-E^xUQ'^)^|«(||g)". (22) 

where the anomalous dimension for the moment is dn = 2 Cpjn/ Po- 

We also show the effect of leading logarithm resummation, (122|) . for a wide range of values 
of Q in Fig. [TJ One can clearly see that this RG-improved result seems to have an even 
minor impact compared with the preceding example ^. Since all the are positive, as 
— i> oo, one expects F — > 0. However, as indicated in Fig. [1], to reach this asymptotic 
result, an unphysically large Higgs mass is required. 

We can also iteratively derive the closed form for the leading logarithms occurring at the 
n-th loop order, by simply replacing the kernel Vq in ( |T6l) by V±. Using ( lA3al) and ( lA3bf) . 
we can predict the leading logarithms through the two-loop order: 



(3o{ ^-21n2j +CpQ-121n2 + 41n2 2 



+ ■ ■ 



As a check, we have explicitly computed the logarithmical contribution in T^^\^) for this 
process It agrees with Cn given above. As one can tell from Fig.[T], even though including 
the leading logarithm at two- loop order noticeably modifies the one- loop result, the overall 
effect of logarithms is too modest to be phenomenologically relevant. 



^ Our investigations on these two examples show that, even when Q 3> m, leading logarithm resummation 
may not play a significant role. This finding provides a counterexample to the assertion made in Q]. 
There is a subtlety arising in this process. Since the composite operator bb acquires an anomalous dimen- 
sion, the renormalization procedure has to be performed to the T^^\ where a renormalization scale iin 
will enter and give rise to the logarithm like ln((5^//i^). This type of logarithm has a very different origin 
from the coUinear logarithm we are interested, so will be eliminated by choosing fin ~ Q. 
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Summary. In this paper we suggest that, for a class of exclusive single-quarkonium 
production processes, the short- distance coefficient in the leading order NRQCD expansion 
can be further separated into "hard" and "coUinear" degrees of freedom. As a consequence, 
the corresponding NRQCD matching coefficient can be expressed as the convolution between 
a hard coefficient function and a perturbatively calculable LCDA of a properly-chosen free 
quark-antiquark pair. The procedure of refactorization has both conceptual and technical 
advantages over the conventional NRQCD matching calculation. For instance, we have 
shown that large kinematic logarithms can be readily summed to all orders in within this 
scheme, which is otherwise difficult to accomplish. This strategy also provides an efficient 
means to reproduce the NRQCD matching coefficients when proceeding beyond the tree 
level (for an explicit one- loop illustration, see [3]). 

We believe that the idea of refactorization is rather general and deserves to be further 
explored in other more complicated cases. For example, the NRQCD short- distance coeffi- 
cients associated with exclusive P-wave quarkonium production and relativistic corrections 
to S-wave quarkonium production, should also be amenable to a similar coUinear factoriza- 
tion formula. It is conceivable that exclusive double charmonium production processes like 
e^e^ J/ip + rjc, may even be approachable from this angle. However, since such processes 
often violate the helicity selection rule and the amplitude usually starts at subleading power 
of m^/Q^, consequently higher-twist LCDAs must be introduced 

Baa 

. This may pose 

some great challenge to fulfill factorization and evolution. 
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Note added. After the manus crip t was submitted, we became aware of a work done by 



Shifman and Vysotsky in 1981 j29j, who had considered two basically identical processes 
as ours, and also attempted to resum the respective leading kinematic logarithms. Though 
differing in technical taste, their and our approaches are essentially the same, both derived 
from the light-cone OPE. Nevertheless, our approach seems easier to follow in practice. For 
the first process e^e~ —>■ ribj, they also considered the mixing effects due to the two-gluon 
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Fock components of rjb, which we have not. Both of works agree on the one-loop coefficient 
of leading logarithms in this process. Our approach enables one to readily deduce the 
two-loop coefficient as well, but it would be a laborious undertaking for their formalism to 
achieve this. For the second process h — > T7, they combined the collinear logarithms with 
those logarithms arising from the running quark mass. After subtracting the latter type 
of logarithms, their results agree with ours. It should be noted, however, the strength of 
refactorization proposed in this work is not only limited to summing large logarithms. As has 
been stressed throughout the paper, our approach will serve an efficient and systematic tool 
to reproduce the NRQCD short- distance coefficients to any desired order in as, including the 
logarithms as well as the constants, in various exclusive quarkonium production processes. 



In this regard, we feel that our work has gone one step farther than Ref. 



29| 



APPENDIX A: USEFUL MATHEMATICAL FORMULAS 



In Example 1, in order to iteratively deduce the leading logarithms through the two-loop 
order from (fT6|l . we need know the following integrals: 



dxT^^\x)Vo{x,y) 







1 [3 

-r, V o + (l-2/)lny + ?/ln(l-?/) 



(Ala) 



dx[ dyT^'\x)Vo{x,y)Vo{y,z) = 

2 ^ , .1„2/ 



9 TT^ 

h\nz + ln(l - z) + 2(1 - z)\nz 

4 6 



+2z\n{l - z) + {l- z) \n^z + zln%l - z) + zLi2(z) 



:i-z)LUl-z 



(Alb) 



where Li2 denotes dilogarithm. Eq. (lAlaO can be found in Ref. 26|], and (lAlbl) is new. 

Expand the Gegenbauer-moment-summation formula ( ]13aP in d2n, truncate the power 
series in as{Q'^) ln(Q^/m^) to the second order, and match them onto Eq. (|T7j) . Enforcing 
the mutual equivalence leads to the following mathematical identities: 



n=0 



1(0) 
'2n 



- -2(3-21n2). 

n=0 

00 2 

,2 1(0) - ^ 



E^2n02n = 9 - y - 16 lu 2 + 2 lu^ 2 . 

n=0 



(A2a) 
(A2b) 
(A2c) 



Analogously, in Example 2, to iteratively acquire the leading logarithms up to two-loop 
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order, we need to know the following integrals: 



dxT^^\x)V±{x,y) 



yii-y) [2 



+ lny + ln(l — y) 



(A3a) 



dx dyT^^\x)V^{x,y)V^{y,z) 

^0 

9 



- + 31n2 + 31n(l-;z) + In^ 2 + ln^(l - ^) 



(A3b) 



z{l-z) 

Expand the Gegenbauer- moment-summation formula (l22l) in d2n, truncate the power 
series in ^^(Q^) ln(Q^/m^) to the second order, and compare them with fl23|) . The mutual 
equivalence demands that the following identities must hold: 



n=0 

oo 



(0) 
n r2n 



•2 ;;(o) 

2n 



-2(3 -4 In 2) , 
9-241n2 + 81n2 2. 



(A4a) 
(A4b) 



n=0 
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